Are Chaplygin gases serious contenders to the dark energy throne? 
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We study the implications on both background and perturbation evolution of introducing a Chap- 
lygin gas component in the universe's ingredients. We perform likelihood analyses using wide- 
ranging, SNla, CMB and large scale structure observations to assess whether such a component 
could be a genuine alternative to a cosmological constant, A. We find that the current data favors 
behavior in an adiabatic Chaplygin Gas that is akin to a cosmological constant. 



m ■ 
o ■ 
o : 

s . 

<: 
oo : 

(N : 
> ■ 

oo : 
o . 

m ■ 

o : 

m . 

o ■ 



Of 
i 

o 

CO 



X 



I. INTRODUCTION 

Supernovae observations 0, Q first indicated that the 
universe's expansion has started to accelerate during re- 
cent cosmological times. This, and further observations, 
e.g. of the Cosmological Microwave Background (CMB) 
or Large Scale Structures (LSS) , suggest that the energy 
density of the universe is dominated by a dark energy 
component, with a negative pressure, driving the accel- 
eration. One of the substantive goals for cosmology, and 
for fundamental physics, is ascertaining the nature of this 
dark energy. Maybe the most attractive option would be 
a cosmological constant, A, however there are infamous 
fine-tuning and coincidence problems associated with ex- 
plaining why A should have today's energy scale. These 
problems have lead to a wealth of dynamical, scalar, dark 
energy ("quintessence") models being proposed as alter- 
natives to A (see Q for a good review). Even in these 
cases, however, explaining why our epoch should be so 
crucial in triggering the acceleration still requires fine- 
tuning. 

A concurrent problem is the nature of the non- 
baryonic, clumping dark matter component required in 
the standard model to give Large Scale Structure predic- 
tions consistent with observations. 

Recently an alternative matter candidate, a General- 
ized Chaplygin Gas (GCG), has been proposed as a po- 
tential 'hybrid' solution to both the dark energy and dark 
matter problems. The GCG can be seen to evolve in a 
wide range of contexts, for example from supersymmetry, 
tachyon cosmologies |4J and brane cosmologies ,5] . A re- 
cent letter |(j dealt with the implications for the matter 
power spectrum in the absence of CDM and effectively 
ruled out the GCG as a CDM substitute. 

In this paper we investigate the strength of the GCG 
as a dark energy candidate. Although there have been 
a number of papers discussing various aspects of GCG 
behavior (0 _ E2) there has not been, as yet, a full anal- 
ysis of the constraints that can be placed on such models 
from the wide range of complementary data sets currently 
available. This is necessary if such exotic matter types 
are to be considered as serious alternatives to the A-CDM 
scenario. 

In section[n]we review the background evolution of the 
GCG and discuss the implications for supernovae (SNla) 



observations. Although such constraints are important, a 
wide range of proposed theories can generate the required 
expansion profile (see [3j for dark energy theories and, for 
example, [l3| for an alternative to dark energy) . In order 
to better discriminate between theories, perturbation de- 
pendent observables must be taken into consideration. In 
section ITTT1 we extend our discussion to perturbations in 
the Chaplygin gas and discuss the implications for struc- 
ture formation in the presence of an adiabatic Chaplygin 
fluid. In section [Tvl we consider the effects on radiation 
perturbations and the CMB spectrum. In section Ivl we 
present the main results of the paper, likelihood analyses 
for a CDM+GCG+baryon universe. We include the op- 
tion of a pure GCG + baryon scenario (fi c = 0) for com- 
pleteness. We obtain a clear indication of the strength 
of the GCG model when compared to CDM and A. In 
section lvTl we summarize our findings and assess the true 
potential of Chaplygin gases as a dark energy contender. 



II. BACKGROUND EVOLUTION 



The Generalized Chaplygin models can be character- 
ized by three parameters: wq, a and f2^. The equation 
of state nowadays w(a = 1) = — \wo\ and the index a 
specify the equation of state evolution, 



\w \n° ch po 



(1) 



where po = 3Hq (8itG = 1) is the total energy density 
today. The energy conservation equation, p + 37Y(1 + 
w)p = 0, admits a solution for p(a) specified by \wq\, a 
and the fractional energy density today, f2° h , 
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The equation of state then evolves as, 



w(a) 
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At early times the GCG's equation of state tends to zero, 
mimicing CDM . The value of a determines the redshift 
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FIG. 1: Contours in fl c h — \wo | space with the same luminosity 
distance as a fiducial model A-CDM model with Qa = 0.7 at 
z = 0.5 (top) and z = 1 (bottom). For a = (full line) the 
luminosity distance curve is identical to that of the fiducial 
model at all redshifts. 



of transition between the two asymptotic behaviors; the 
greater the value of a the lower the transition redshift. 
At early times, the total amount of matter with w ~ 
reaches an asymptotic value 
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(4) 



where Vl m is the baryonic + CDM density fraction. Note 
that the unique ability of the GCG to account for both 
the dark energy like behavior at late times and for ordi- 
nary dark matter at early times motivated the original 
studies of this particular equation of state. 

In previous discussions a has often been assigned a 
positive value in the range < a < 1, in order to be 
consistent with various higher dimensional theories that 
can produce a perfect fluid stress-energy tensor satisfy- 
ing the criterion in £[J (see for example |14|). In our 
study, we extend the range of values of a considered to 
— 1 < a < oo in order to obtain a broader assessment of 
whether Chaplygin gases could be a viable alternative to 
the standard model. 

For a=0 the background evolution of the Chaply- 
gin gas is identical to a A-CDM model with fiA,e// = 
Q c h\w \, andfi m) e// = n m + n c fc(l- \w \)- Furthermore, 
as is visible in equations an d ©, when |tuo| tends to 
1, the GCG component tends to evolve as a cosmological 
constant, irrespective of the value of a. Note that there is 
no analogous quintessence like behavior (with wq —1), 
thus we are only comparing GCG to theories including 
A. 



The SNla observations measure the apparent magni- 
tude, m(z), related to the luminosity distance, cLl{z) via, 



m(z) = M + 5\ogd L (z) + 25, 
dz' 



d L = (l + z) f 
Jo 



H(z') 



(5) 
(6) 



where M. is the absolute bolumetric magnitude and 
is measured in Mpc. It is easy to see that, because the 
background evolution (through H) wholly determines lu- 
minosity distance predictions, the degeneracy between a 
GCG with a = and A-CDM will allow the Chaplygin 
gas to fit the SNla data well. Indeed the degeneracy also 
stretches to a ^ when one considers luminosity distance 
at a specific redshift. In figure^we show Chaplygin mod- 
els with degenerate luminosity distances with a fiducial 
A-CDM model with O A =0.7, at z = 0.5 and 1.0. This 
degeneracy, however, implies that the SNla observations 
cannot be a strong discriminant between the GCG and 
A; we must look to alternative, perturbation-dependent 
observations to test the validity of the GCG models. 



III. CHAPLYGIN GAS PERTURBATIONS 

We treat the Chaplygin gas as a perfect fluid made 
up of effectively massless particles interacting with the 
rest of matter purely through gravity. We assume purely 
adiabatic contributions to the perturbations so that the 
speed of sound for the fluid is 
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(7) 



and the time variation of w is 



w = -m (1 + w) (c 2 s -w)= 3Hw(l + w)(a + 1) (8) 

where derivatives are with respect to conformal time 
(d/dr), and aH = da/dr. 

In the synchronous gauge and following the approach 
and notations of Ma and Bertschinger |15| , we can write 
down the evolution equations for the density and velocity 
divergence perturbations, S and 9, using the conservation 
of energy momentum tensor T^.^ = 0, 



5 = -(1 + w) \6+- \ -3H(c 2 s -w)5, (9) 



= -W(1-3<S). 



(l+w) 



k 2 S - k 2 a. (10) 



The fluid is highly non-relativistic and therefore we as- 
sume the shear perturbation a = 0. 

At early times, when the Chaplygin gas has w w 0, 
the GCG perturbations evolve like those of ordinary dust 
with 9 = 8 = 0, and S = —h/2. In the radiation era 
5(a) oc a 2 , while S oc a in the early GCG dominated era. 
At later times, when the GCG's equation of state starts 



3 











: ' ' k/H-'ioo ' : 


k/H='lOOO ' ' : 


\ — i — i — i — I — — i — I — i — i — - 


i — i — i — i — i — — i — i — i — i — - 


: k/H=i : 
: 1 [3] : 

: [1] : 


k/H=10 : 







0.5 It 0.5 



l w ol 

FIG. 2: Late-time 5 evolution as a function of a and |wo| 
for 4 scales k/H = 1,10,100,1000 [H « 3000ftA'/pc _1 , h = 
0.65). Regions [1] and [2] undergo power law growth and 
decay respectively and [3] undergoes oscillatory decay. 
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FIG. 3: Late time evolution envelope for S. The power, 
7i e // = tS/8, is plotted for two length scales k/H — 10 and 
100. In the oscillatory regime n e / / for the bounding envelope 
is plotted. 



to decrease, the perturbations stray drastically from this 
dust-like evolution. 

We can understand the late-time behavior more clearly 
if we evaluate the second order differential equation for 
6. By differentiating equation 10 with respect to time 
we find, as outlined in an appendix (section IVIIj) , for a 
general, shear less, fluid, 

5 + [1 + 6{c 2 s - w)] HS 



9H 2 (c 2 s ~ wf + 3H(c 2 -w) + 3-( C 2 - w) + c 2 s k 2 



-3c 2 (l + w)W + — (1 + w)(38P + Sp). 



(11) 



Numerical integration shows that the coupling to 9 in 
equation (|ll|l is subdominant for all scales that we are 
interested. For the Chaplygin gas, 

5 + AHS + Bn 2 5-^(l + w)n c S c = 0, (12) 



A = 1 - 6w(a + 1), 



(13) 



{fl ch + [7 + n ch + (13 - 3Q c/ > + 6a 2 
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-3n ch (l + 2a)w 2 
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(14) 



where the subscript c refers to cold dark matter. 

For w = 0, ijT2|) reduces to the expected, scale indepen- 
dent, CDM evolution with 5 oc r 2 . For io^0we retain 
scale independence if a = and just get suppression of 



density perturbations. However for a/0 the perturba- 
tion evolution becomes scale dependent with the (k/H) 2 
term dominating the others for scales greater than a char- 
acteristic scale 



u 2 



\aw\ 



There are 3 possible solution types, for k 3> A;*, 



[1] 
[2] 
[3] 



growing mode, a < 0, 
decaying mode, a ~ 0, 
oscillatory decay, a > 0. 



(15) 



(16) 



In figure[21we show the asymptotic behavior at late times 
(taking w ~ — \wq\) as one increases the scale k/H. In 
figure |3| we show the associated scaling of 5, plotting 
n eff = tS/5. For \w \ ^ 0, and a > 0(< 0) GCG per- 
turbations are suppressed (promoted) in comparison to 
those for a A-CDM model. 

The GCG also has an effect on the the CDM pertur- 
bations through the relation: 



XH 2 

H5 C - [n c 6 c 



+ (1 - 3aw)Ct ch S] = 0. (17) 



If \w \ 7^ and a > 0(< 0) the GCG drives suppression 
(growth) in S c . In hgureQ]we show the power law evo- 
lution of S c , plotting n e ff yC = tS c /S c for 3 scales as one 
varies a. The strong growth in S c effectively rules out a 
GCG with a < as a dark energy candidate. 

Note that the a=0 degeneracy present in the back- 
ground evolution is not found in the perturbations. The 
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FIG. 4: Evolution of n e //, c = tS c /5 c for \wo\ = 0.5 and a=- 
0.1 (short dash), (long dash), 0.1 (long-short dash), 0.5 (dot 
dash) in comparison to pure CDM model (full line). Three 
length scales k = 10~ 4 , 10" 2 , W^hMpc' 1 are considered. 



FIG. 5: Large scale ISW effect for GCG with q=-0.5 (long 
dash), (short dash-dot) and 1 (long dash -dot) compared 
against a CDM+baryon (short dash) and 'fiducial' S7a = 
0.7, Sl m = 0.3 (full line) models. 



matter power spectrum and CMB spectrum will there- 
fore be better discriminators between A and the GCG 
than SNla. 



IV. IMPLICATIONS FOR TEMPERATURE 
ANISOTROPIES 

A Chaplygin gas matter component would change the 
temperature anisotropy specturm in a number of ways; 
altering the late-time ISW effect, the peak positions and 
relative heights. 

The GCG's late-time evolution will alter the evolu- 
tion of the gravitational potential the CMB photons pass 
through to reach us, inducing an ISW effect. Following 
[l6| and again using the terminology of [T^. the ISW 
temperature anisotropy is given by a source, 



S 



ISW 



(18) 



dr 



--{p + P)a 2 (j - a 2 Sp - 3Ha 2 (p + P)-- 



where $ and $ are the Bardeen variables [l7| . 

At late times the shear, a, is negligible and it is the 
density perturbation that drives the ISW effect. 



I [-"*] 



(1 + 3™) - 



l eff 



Hp5 (19) 



where n e ft is the power law index and p — fir as de- 
scribed in section IIIII Equation i|19|) shows why in a 
standard CDM scenario, with w = and n e ff and p 



both ss 2, there is no appreciable ISW effect. In section 
IIIII however we saw for a < that n e f f > 2 giving a 
negative ISW effect, while for a > 0, n e // < 2 producing 
an increase in the ISW temperature anisotropy. These 
effects are shown in figure 

The position of the first peak will be altered through 
adjustments to the sound horizon, r s h or , and angular di- 
ameter distance at the last scattering surface, cIa- The 
position of the first peak in multipole space is given by 



u = 



nd A {z rec ) 

Tshor (Zrec) 



with 



dA 



H 



(20) 



(21) 



da 



Tsho 



cfdT 



(22) 



m.eff 



a — cl b da 



a 2 



where e ^ is defined in equation |@J and c] b is the 
speed of sound for the radiation-baryon system, not to 
be confused with c s for the Chaplygin gas (at this time 
the Chaplygin gas is behaving like dust and has c 2 = 0). 
For fixed u>b = flbh 2 and uj c = Q c h 2 (h = H o /100), I a 
increases as one increases a or \wq\. The position of the 
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FIG. 6: Comparison of CMB power spectra, normalized to 
COBE at I = 10, with varying a and |wo|, keeping all other 
relevant quantities fixed. The behavior of peak heights and 
positions is discussed in the text. 



first peak would be the same for scenarios with the same 
value of ^l i eff d A . 

Because the Chaplygin gas mimics matter at early 
times there is no simple degeneracy, governed by the 
peak positions, as there is for quintessence models (see 
for example (HI)- The peak heights, when compared 
to the low I 'plateau', depend upon a and \wq\ through 
their influence on the ISW effect, the horizon scale at 
matter-radiation equality, l eq , (through Cl^^fh ) and 
the depth of the potential well at last scattering (also 
through n° rneff h 2 ). 

We follow the phenomenological discussion in [19|, to 
predict how the peak heights will alter for fixed ujf, and u> c . 
Increasing a, increases ^„ e j^, so that matter-radiation 
equality happens earlier, increasing l eq and curtailing the 
driving effect that the decay of the gravitational poten- 
tial has on (5 7 oscillations during the radiation era. This 
lowers the height of the first peak, a decrease which is 
compounded by the raising of the plateau from the ISW 
effect. An earlier matter-radiation equality also decreases 
the depth of the potential well at last scattering, which 
combined with the reduction in radiation driving, in- 
creases the height of the third peak in comparison to 
the first and second ones. 

As one increases \wq\ one decreases Cl^ ne fp lowering 
l eq , and increasing the height of the first peak. This 
is tempered, however, by the increase in plateau height 
from the ISW effect. Reducing l eq acts to decrease the 
height of the third peak in comparison to the second and 
first ones. These behaviors are confirmed by the full anal- 
ysis, as is shown in figure |SJ 



The multitudinous effects that the GCG has on the 
CMB spectrum make comparison with CMB observa- 
tions a strong test for the GCG models as will be seen 
below. 



V. CHAPLYGIN GAS LIKELIHOOD ANALYSIS 

In order now to assess the viability of a 
GCG+CDM+baryon universe, we turn to evaluate 
the probability (the posterior) of these models given 
some current observations, namely SNla, CMB and LSS 
probed through galaxy survey. 

To study the posterior distribution, we use the Baye's 
theorem and rewrite it as the product of the likelihood 
and the prior (we assume the evidence is constant and 
thus ignore it). To probe this posterior, we consequently 
compute both the likelihood and the prior at various po- 
sitions in the chosen restricted parameter space. This 
sampling is conducted via the construction of a Monte 
Carlo Markov Chain through the Metropolis-Hasting al- 
gorithm. Once converged, this chain provides us with a 
collection of independant samples from the posterior (see 
[20l I2H 122^ for an introduction to this technique in this 
context and (2^, 1^3 f° r general guidance). 

Our code uses some likelihood computation elements 
from the code described in [22J , and relies on a version of 
the CAMB code [25| extended to include a Chaplygin gas 
component in order to calculate CMB power spectra and 
matter power spectra. As input data, we considered the 
apparent magnitudes of 51 Supernovae @, CMB data 
sets from COB E El , MAXIMA [27|, BOOMERANG 
|28j . and VSA |29| and large scale structure data from 
2dF [3(j. We consider only flat models, i.e. VLk = with 
scale invariant initial power spectrum, i.e. n s = 1. We 
use stringent (Gaussian) priors on Hq using the HST Key 
Project results h = 0.72 ± 0.08 UK and on uj b = Q b h 2 = 
0.02 ± 0.001 using BBN constraints pi] . 

We normalize the matter power spectrum using A s , the 
initial power spectrum normalization, and following |33j , 
we use and b± to parameterize redshift-space distor- 
tions and (linear) bias respectively. The power spectrum 
is then related to the transfer function T(k) (computed 
with CAMB) by 

^)=^(l + y + y)^T(fc) 2 . (23) 

In order to alleviate the natural degeneracy between A s 
and bi (as far as LSS constraints are concerned), we use 
the 2dF results 0,113 to impose strong (Gaussian) pri- 
ors on (3 and 61, i.e. f3 = 0.54 ± 0.09, b x = 1.04 ± 0.11. 

Throughout this analysis, we ensured the chains' con- 
vergence by generating and comparing several of them 
(typically containing 10 5 elements) and by checking the 
so-called "parameter mixing" amid them. After several 
trials, we choose the proposal density for each parameter 
to be a Gaussian whose width is close to the final one 
and whose center is the last chain values. This allows a 
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FIG. 7: Joint posterior of the a and — Wo parameters considering only SNla data (top- left panel), LSS data (top-right panel), 
CMB data (bottom- left) and jointly CMB and LSS data (bottom-right). The contours represent the subsequent 68% and 95% 
confidence regions. While SNla data induce constraints that are quite loose, LSS and CMB constraints are much tighter and 
tend to favor a cosmological constant like scenario. 



full exploration of the parameter space. To pick-up the 
next chain element, we allow only 1 to 3 directions (this 
number is randomly chosen) to vary. This gives us an 
acceptance rate around 25%, a good target value for ef- 
ficiency's sake [i^. The first 4000 elements of the chain, 
prior to its convergence, are thrown away and no extra 
thinning is applied 23] . 

Once converged, the chains provide a fair sampling of 
the full posterior distribution so that we can deduce eas- 
ily from it all the quantities of interest, e.g. the (joint) 
marginalized distribution of any parameter (s). 

As stated above, we are interested in finding the com- 



patibility of a Chaplygin gas + CDM + baryon universe 
with current data. For this we vary only 8 parameters 
{h, Lib, LOcdm, ot, wo, b, (3, A s } and impose the priors stated 
above. We allow a free proposal distribution for u) c , in- 
cluding lo c = consistent with a unified matter universe 
purely containing a GCG and baryons. This allows the 
full breadth of GCG roles (as both a dark matter and 
dark energy candidate) to be tested. Following section 
IIIIl s discussion, we restrict ourselves to < a < 1 and 
— 1 < wq < 0. In figure we plot the marginalized 
joint distribution of the a, \wq\ parameters (which is in 
this case just the joint number density of those param- 
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eters), as well as the 68% and 95% confidence contours, 
considering separately SNla, LSS, CMB data sets, and 
also jointly CMB and LSS. Note that for visual purposes 
only the displayed surface has been build by oversam- 
pling our samples using cubic interpolation. This does 
not affect the quantitative interpretation since the distri- 
butions turn out to be smooth. 

The interpretation of the contours is nicely consistent 
with the theoretical prospects discussed above. First, the 
SNla observations (top-left panel) offer very light con- 
straints on the GCG parameters, since they are sensitive 
only to the background evolution. Any a value appear 
viable, extending thus the obvious degeneracy between 
A-CDM model and GCG models with a — discussed in 
sectionllTlf e.g. see figure^. As soon as density perturba- 
tions are considered, the constraints tighten drastically. 
For both LSS and CMB, the isocontours are roughly cen- 
tered on the a = 0, wq = — 1 model, that corresponds to 
the GCG acting like a A term. This fact is emphasized 
in the joint CMB + LSS analysis. Note that in the limit 
that w tends to —1 (we however impose |u> | < 1), the 
GCG component tends to behave like a A term, irrespec- 
tive of the precise value of a, thus leading to the observed 
degeneracy in the — wq — 1 direction. 

The other varied parameters, i.e. h, u>b and Lu c dm, as 
well as the flatness imposed Cl c h, exhibit (joint) distri- 
butions similar to those found in typical A-CDM model 
studies (see e.g. [l^])- This leads us to the the main con- 
clusion of this study: the current data tends to favor ordi- 
nary A theory. When marginalized over all other parame- 
ters, we indeed find, a < 0.5, 0.93, and w < —0.85, —0.8, 
both respectively at the 68% and 95% confidence level. 

VI. CONCLUSIONS 

We have investigated the effect of a Chaplygin gas mat- 
ter component in the universe's ingredients, to see if such 
a component is consistent with observations and whether 
it is a feasible alternative to CDM and A. 

Through inherent degeneracies with A in the back- 
ground evolution, the Chaplygin gas models have a good 
fit with SNla data. These degeneracies are not present, 
however, in the perturbation evolution. In particular 



the growth/suppression of both GCG and CDM density 
perturbations proves distinctive when comparing against 
large scale structure observations; this statement is valid, 
of course, for all cases (for all a) except, naturally, 
wq— = —1 which is identical to A, and has no per- 
turbations. The GCG also introduces a number of dis- 
tinguishing differences from the A-CDM CMB spectrum 
through altering the potential at last scattering, the ISW 
signature, the equality scale, and the angular diameter 
distance to last scattering. Combined, these differences 
provide a strong test for the GCG scenario. 

We performed likelihood analyses using SNla, CMB 
and LSS datasets and found that the current data 
strongly prefers a A- like dark energy component, with 
a < 0.5 and w < -0.85 at the 68% level and with CDM 
as the preferred pressureless matter component. Note 
that, in comparison, the unified dark matter model, with 
fl c = 0, is highly disfavored by the data. This result 
is consistent, but considerably tightens, previous con- 
straints from supernovae, CMB peak position and matter 
power spectrum shape parameter analyses (0 _ E3)- Our 
constraints can be recast in terms of the 'statefinder' pa- 
rameters of |H, r < 1.20 and s > -0.075 at the 68% 
level, thus greatly reducing the ability of a Chaplygin 
Gas to explain the 'cosmic conundrum' problem as pro- 
posed in |37| . 

Our analysis assumed adiabatic perturbations for the 
Chaplygin gas; it remains to be seen how enriching this 
model by considering non-adiabatic perturbations, as 
mentioned in a paper presented after the initial posting 
of this work [3^], might alter the analysis. 

On the basis of current observations however, Chaply- 
gin gases, with adiabatic perturbations at least, do not 
seem to provide a favored alternative to scenarios involv- 
ing CDM and a cosmological constant. 
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and equation of state and speed of sound equations: 
5P P 

c 2 = — = — (assuming adiabaticity) (25) 
dp p 

= -3H{l + w)(c 2 s -w) (26) 

Following [15j the first order perturbation equations 
are 



5 = -(1 + w) U+- \-3(c 2 s -w)m (27) 




(l-3c 2 ) W 



(1 + w) 



k 2 S. 



(28) 



So that the second order equation in 5 (differentiating 
© is given by 



S = -{1 + w) 

-3 (c 2 s - w) (H6 + H5) -3[c 2 -w) H5. (29) 




We eliminate the time derivatives of the metric pertur- 
bations, h and rj, using the perturbed Einstein equations 



k rj 
h 



1 



-Hh 



1 



-a Sp 



2H - 2k 2 rj = -3a 2 SP 



(30) 
(31) 



which give 



(1 + w)- = HS+(l + w)H0 + 3(c 2 s -w)H 2 5 
a 2 

-(l + w) — (S P + 3SP). (32) 

Collecting terms together we obtain the general evolu- 
tion equation for 5 for any fluid with equation of state w 
and speed of sound c s , 

<5 = -3c 2 s (l + w)U6 - [1 + 6{c 2 - w)]HS - [c 2 k 2 + 



9(cf - w) 2 H 2 + 3[c 2 s -w)H + 3-(c 2 s -w) 



VII. APPENDIX - PERTURBATION 
EVOLUTION FOR A GENERAL, ADIABATIC 
FLUID 



We use the basic background equation (time deriva- 
tives with respect to r): 



u + n 2 ='- = ^ 



1-3^(0^) 



(24) 



+ (l + w) — (5p + 35P). 



(33) 



Specializing to the Chaplygin gas in the matter domi- 
nated era 



-(6p + 36P) 



3H 2 

2 
H 2 



[n ch (l + 3c 2 s )5 + n c S c ] (34) 



- = -^-(1 - 3fl ch w) 
a I 

„2 



(35) 

aw (36) 
3aHw(l + w)(l +a) (37) 



9 



we find 

<5 + [1 - 6w(a + l)]H5 - 



awk 2 + 



m 2 



Similarly if one applies equation l|33|) to CDM with 
w = c 2 = 0, we find 



{n ch + (7 + n 



ch 



+ (13 - 3Q ch )a + 6a 2 ) w - 3Q c/[ (l + 2a)w 2 }] 6 
3H 2 

—-(l + w)tt c 5 c + 3aw(l + w)H8. (38) 



I + ns c 



m 2 



[Q C S C + (1 - 3aw)Ct ch 6] = 0. (39) 



